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Abstract 

Hadron production at forward rapidities in d+Au collisions is studied in the frame- 
work of parton recombination in the final state. Multiple scattering and gluon satu- 
ration in the initial state are not explicitly considered. The recombination of soft and 
shower partons is found to be important. The soft parton densities are determined 
by extrapolation from the parametrization fixed at r/ = with no unconstrained ad- 
justments. The suppression of Rcp observed at high r\ is understood as the simple 
consequence of the reduction of the soft parton density on the deuteron side compared 
to that on the gold side. The effect of momentum degradation responsible for baryon 
stopping is also considered. The asymmetry of backward-to-forward spectra can be 
reproduced within the same framework without any change in the basic physics. 
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1 Introduction 



Particle production at forward rapidities in d+Au collisions has recently been recognized 
as the fertile ground for testing models of hadron production that embody diverse physical 
mechanisms. Since no hot and dense medium is created in d+Au collisions, they have gener- 
ally been regarded as the type of collision processes where the effects of final-state interaction 
(FSI) are minimal, thus allowing the physics of initial-state interaction (ISI) to manifest itself 
most transparently P0-[H|. That is especially true at forward rapidities where the saturation 
effect is expected to be important and the physics of color glass condensate (CGC) may 
dominate PP-jHj- Alternative approaches to the problem based primarily on perturbative 
QCD (pQCD) have also been considered with emphasis on the effects of nuclear shadowing 
pHIH|. The results do not provide unambiguous interpretation of the data from BRAHMS 9J. 
We consider in this paper yet another approach that emphasizes the hadronization part of 
the FSI, and takes into account hard scattering in pQCD; however, instead of fragmentation 
we consider the recombination of soft and shower partons. Our treatment is an extension 
to the forward region of what has been found to be successful at midrapidity for Au+Au 
collisions ^Hj and for d+Au collisions [TT] . and is a natural mechanism for hadronization at 
intermediate transverese momentum p?- It does not, in principle, contradict the physics of 
CGC. The issues are whether hadronization by recombination is important and which one of 
the competing effects is dominant in the kinematical region under examination. Our results 
indicate that the BRAHMS data jHj can be well reproduced in our approach even when no 
adjustable parameters are used. 

Experiments at the Relativistic Heavy Ion Collider (RHIC) have found that the ratio, 
Rcp, of d+Au central to peripheral inclusive spectra for 1 < p T < 3 GeV/c decreases 
monotonically from a value ~ 1.8 at pseudorapidity rj ~ — 2 to a value ~ 0.5 at 77 ~ 3.2 
H21- This has led to the interpretation of a change of the physics responsible for the 
phenomena from the gold side (rj < 0) to the deuteron side (rj > 0) [I]. For 77 < the 
enhancement of the particle yield at intermediate px with respect to binary collision scaling, 
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i.e., the Cronin effect [T3], is generally regarded as the result of multiple scattering in ISI 
[THUS]. For V > saturation physics is considered to be dominant, especially at large 77, so 
that there is suppression, instead of enhancement, in particle production 00111. Neither 
effect takes into account any details about hadronization in FSI. Such explanations of the 
data cannot account for the phenomenological fact that R V CP (for protons) is larger than 
R^p (for pions) at intermediate 16]. Hadronization of partons to a proton or a pion is an 
issue that involves the final state, and the use of fragmentation functions (FF) appropriate 
for the produced hadrons in the usual factorizable way would necessarily yield a p/ir ratio 
< 1 by virtue of the nature of the FFs, contrary to what has been observed [TT] . 

The subject of large p/n ratio has been addressed by several groups that consider parton 
recombination as the hadronization process in the final state in Au+Au collisions [T%1 IT9"1 12U] . 
For hadrons produced at intermediate pt it is the recombination of thermal partons at lower 
Pt that gives rise to the higher yield of the baryons compared to mesons. In Ref. [113] the 
additional component of the recombination of thermal and shower parton is introduced, and 
is shown to be dominant in the 3 < pt < 9 GeV/c range. That component turns out to 
be crucial to explain the Cronin effect, when the formalism is extended to d+Au collisions, 
without the need for any fey broadening in ISI jllj . Although there are no thermal partons 
in d+Au collisions as in Au+Au collisions, there are soft partons that play the same role. It 
is found that the property Rqp > Rcp a ^ intermediate px can readily be reproduced in the 
recombination model fUE]. Based on the success of that description of particle production 
at midrapidity in d+Au collisions, we now extend the treatment to the forward region (77 > 0) 
without the introduction of any new physics. This extension should be considered whether 
or not the signature of new physics is present in the forward region, since the contribution 
from conventional physics forms the background that must be understood to facilitate the 
identification of any new signal. As we shall show, it seems that recombination in FSI is 
completely adequate to explain the suppression at large 77. 

The main reason for Rqp to decrease with increasing rj is simple. Since the recombination 
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process that involves soft partons is important, the dependencies of the soft parton density on 
rapidity and centrality directly affect the particle spectra at intermediate px- That density 
that is reflected in dN/ drj is known to decrease with increasing rj with the consequence that 
Rep also decreases with 77. Since CGC has successfully reproduced the hadron multiplicity 
dN/dr], one may take the view that the suppression at large 77 has its origin in saturation 
physics. That may be the case. From the point of view of FSI the suppression is a direct 
result of parton recombination involving soft partons whose density decreases with rj; indeed, 
our treatment of hadron production is compatible with any model that can generate the soft 
parton spectra that we shall determine phenomenologically. 

Before we consider hadron production quantitatively, let us make some general remarks 
here on recombination. Since the formation of a pion by a qq pair implies the loss of degrees 
of freedom, it is sometimes thought that the recombination process leads to a decrease of 
entropy. However, the entropy principle should not be applied locally. A global consideration 
must recognize that the bulk volume is increasing during the hadronization process, and 
thus compensates the decrease of entropy density that is a local quantity. Furthermore, in 
the qq — > 7r process color singlet is formed by color mutation of the quarks through soft 
gluon radiation that carries the extra degrees of freedom without significantly altering the 
local relation that the sum of the parton momenta equals the pion momentum. All gluons 
hadronize by conversion to qq pairs first, so there is a cascading of the degrees of freedom to 
lower and lower momentum region where the pion multiplicity accumulates. Eventually all 
degrees of freedom are converted from partons to hadrons. Such non-perturbative processes 
cannot be calculated, and our formalism does not contain explicitly the feed-down process 
of partons cascading to the very low pt region before recombination. Indeed, such a process 
need not be made explicit in our approach, since we do not determine the soft parton 
distribution in a model of evolution from the initial state. Instead, we determine the soft 
parton distribution from the observed pion distribution in the 0.5 < pr < 2 GeV/c region. 
The soft partons thus obtained are defined in the context of recombination. We use them 



4 



in the same context when we consider their recombination with shower partons. Thus our 
procedure of treating hadronization at all > 0.5 GeV/c is totally self-consistent. 

It is also appropriate at this point to make another general remark about our treatment 
of recombination. In a one-dimensional (ID) formulation of the recombination process in 
momentum space, one may question whether we are ignoring the spatial extent of the re- 
combining subsystem normal to the collinear momentum vectors of the quarks and hadrons. 
Partons that are separated by a distance of the order of the transverse size of the bulk volume 
but have parallel momenta are not likely to recombine. That is indeed an issue that has to 
be faced in a model in which the spatial and momentum distributions of the partons are 
generated from the initial state and evolve according to some sensible dynamics. However, 
that is not what we do, as we have already stated at the end of the preceding paragraph. For 
the soft parton distribution, which is the crux of the issue here, we start from the observed 
soft pion distribution. For every p of such a pion we claim that it can only arise from the 
recombination of a quark at p\ and an antiquark at p 2 that are not only collinear in mo- 
mentum vectors, but also spatially overlapping within a transverse space of the order of a 
hadron. In other words, since we determine the soft parton distribution from the final state, 
the result automatically implies that only those partons can recombine to give the observed 
hadron. The ID formulation of recombination is therefore appropriate for the way in which 
it is applied. 



2 Pion Production at 77 = 

We summarize first the production of pions by parton recombination at midrapidity in d+Au 
collisions and make minor adjustments to align our calculation for the centrality cuts of the 
BRAHMS experiment [Oj. We use the formalism for hadronization described in Refs. [TTHlllj. 
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where the inclusive distribution of pion at r] = in a ID description is given by 

dN n r dp 1 dp 2 

P—^—= / Fqq>{Pl,P2)Rx{Pl,P2,P), (1) 

dp J pi P2 

p being in the direction of the detected pion. F q qi(pi,p2) is the joint distribution of a q and 
a q' at p± and P2, and R n (pi,p2,p) is the recombination function for forming a pion at p: 
Rtt(pi,P2,p) = (PiP2/p)$(Pi + P2 —p)- For p in the transverse plane so that p T = p, the 
distribution dN n /d 2 pdri\ v=0 , averaged over all (ft, is 



dN n 1 fP 
pdp p 3 



rP 

/ dpiFtffaiP-pi). (2) 
JO 

F q q has three components 

F ?r = rr + rs + ss, (3) 

where T stands for soft parton distribution and S for shower parton distribution. For Au+Au 
collisions T would refer to the thermal partons ^0], but in d+Au collisions the notion of 
thermalization is inappropriate, but soft partons nevertheless exist. Since they are treated 
in the same way, the same notation is used with T, which can be regarded as referring to the 
last letter of "soft." At low pt the observed pion distribution is exponential, which suggests 
the form 

T{Pi) = Pi-r~ = C Pi ex v(-Pi/T), (4) 
dpi 

so that the TT component in Eq. (JHJ) yields 

dN? r _ 
pdp 6 



exp(-p/T) . (5) 
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The values of C and T have been determined in Ref. jTTj already from the low-p-r d+Au data. 
Since our approach to hadron production at intermediate pt at RHIC is to emphasize parton 
recombination in the final state, we use the phenomenological input for the soft component, 
without relying on any specific model for soft partons so that our result will be independent 
of the reliability of such models. 

The distribution S is a convolution of the hard parton distribution fi(k) with transverse 
momentum k and the shower parton distribution (SPD) Sj (z) from hard parton % to semi- 
hard parton j 



where k min is set at 3 GeV/c, below which the pQCD derivation of fi(k) is invalid. For 
each i, fi(k) depends on the parton distribution functions, nuclear shadowing, and hard 
scattering cross sections. The result is presented in the power-law form, whose parameters 
are tabulated in Ref. [TT]. The SPDs are obtained from the FFs and are given in Ref. [22]. 

It is sufficient with the above specification of T and S to calculate the p? spectrum of 
pion by use of Eq. (J2J). However, for comparison with the BRAHMS data we need the values 
of C for the corresponding centrality cuts. Let us use (3 to denote the centrality cut, which 
is an experimental quantity related to the impact parameter h. In Ref. [11J it is found that 
C(p,rj) varies with (3 at rj = according as C((3, 0) = 12,11,7.8,5.65 (GeV/cT 1 for (3 = 
0-20, 20-40, 40-60, 60-90%, respectively. The centrality cuts of BRAHMS |] are (3=0-20, 
30-50, and 60-80%. We therefore make the interpolation and set C((3,0) = 12,9.0, and 6.55 
(GeV/c) -1 , respectively. The value of the inverse slope is T = 0.208 GeV/c, as determined 
in Ref. |TT] . 

Before showing the result, we remark that the SS component in Eq. (j3J) corresponds to 




(6) 
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fragmentation, if the two shower partons originate from the same hard parton, i.e., 

= EjL**M*> (j^)} • m 

where the curly brackets signify the symmetrization of the leading parton momentum fraction 
I22j. If the two shower partons are from two independent jets, then there would be two 
fi(k) distributions. The recombination of such partons is very unlikely to occur in d+Au 
collisions at 200 GeV/c and will be ignored. 

The results of our calculation for tt + production are shown in Fig. 1 for the three centrali- 
zes, displaced by factors of 10 2 from neighboring ones. The light solid lines show the soft-soft 
(TT) components that are straight lines in the log plot. The dashed lines show the soft- 
shower (TS) contributions and the dash-dotted lines the shower-shower (SS) contributions. 
The heavy solid lines are the sums, whose deviations from the straight lines are indicative 
of the effects of hard scattering. Note that the TS contribution becomes less important as 
(3 increases because C(f3,0) decreases. Indeed, if the /3=60-80% case is regarded as being 
almost like the pp collision, we see from Fig. 1 that the neglect of the TS contribution does 
not constitute a bad approximation, and the large pt behavior is essentially governed by jet 
fragmentation, as has traditionally been used to treat pp collisions. However, that is not the 
case for central d+Au collisions. The TS contribution to the spectra is what accounts for 
the Cronin enhancement at intermediate pt without kx broadening due to ISI [TJ] • The pt 
distributions have been shown to agree with the PHENIX data JH]- The Rcp ratio for the 
(3 values of the BRAHMS data will be shown in Fig. 5 below. 

3 Pion Production at 77 > 

Since the formalism for pion production at r] — described in the preceding section success- 
fully reproduces the experimental pt spectra at all centralities, we now make a straightfor- 
ward extension to the 77 > region. This extension should be made with no change in the 
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Table 1: Values of C{f3,rj) in (GeV/c)" 







V 






1 2.2 3.2 




0-20% 
30-50% 
60-80% 


12.0 11.1 9.01 7.05 
9.0 8.5 7.9 6.0 
6.55 6.6 6.1 5.1 



basic physics underlying the formalism in order to provide a baseline for comparison with 
the data before the search for the signature of any other physical origin. 

The quantities that must be modified for the rj > region are C(/3,r)) for T and fi(k) 
for S. For the inverse slope T we shall proceed in two steps. First, we keep T fixed as 
r) is increased, since no data on low-p^ 7r + spectra are available to serve as our guide for 
its modification. It is of interest to see how close the calculated Rep will turn out to be 
in comparison to the data, when the constant-T assumption is applied for the purpose of 
introducing no adjustable parameter. Later, we shall allow T to depend weakly on (3 and rj 
and show that the fit of Rcp can be improved. 

Since the observed rapidity density dN c ^/drj is an integral over the px distribution, which 
is dominated by the soft contribution at low pt, we see from Eq. (J5J) that dN^/dr] should 
be proportional to C 2 (f3,r)). We can therefore determine C(j3,rj) by use of the formula 



C(f3, V )=C(f3,0) 



dN ch /d V ((3) 11 2 
dN ch /dr]\ v=0 ((3) 



C(/3,0) is given in the preceding section, while dN c h/dr)(f3) is known from PHOBOS data 
We thus obtain the values of C(/3,T]) as shown in Table I. Since PHOBOS does not 
have exactly the same centrality cuts as BRAHMS, some interpolation between neighboring 
values have been made to deduce the numbers in Table I. 

We have parametrized the hard parton distributions fi(kr,y), as before [HJI21I, m the 
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power-law form 

for y < 1. For larger values of y the spectra are increasingly suppressed at high kr because of 
the phase space boundary that requires kx < ko(y)- The kinematic limit is given by ko(rj) = 
\/s/(2 coshy) and takes the values 21.9 and 8.13 GeV/c for y = 2.2 and 3.2 respectively. 
To take into account the change of fi(kT,y) from positive to negative curvature around 
k T ps 0.5 k , equation Q must be modified by a dampening factor. Therefore we use the 
parametrization 

(7 \ -rij / , \ mj 

±1 + ^) l 1 -^) > ( 10 ) 

for y > 1 where the ±1 sign is used in accordance to whether the given value of Bi is preceded 
by a + or — sign. 

The parametrizations are obtained from leading order minijet calculations using CTEQ5 
parton distributions [22]. EKS98 shadowing was used for the Au nucleus while the 
deuteron was treated as a superposition of a proton and a neutron without further nuclear 
modifications. The values of all the parameters for central d+Au collisions (J3 = 0-20%) are 
given in Table II, corresponding to a m = 40.3 mb and the K factor unspecified. For other 
centralities scaling in the number of binary collisions, iV co ii, is assumed. We set K = 2 in 
our calculation below. 

Using C(P, rj) in Eq. (JH) and /i(/cr, y) in Eq. ©, neglecting the difference between r\ and 
y for rj > 1, and fixing T at 0.208 GeV/c as for rj—0, we can now calculate dN n +/pdpdr] 
according to Eqs. (0) and © and obtain the results shown in Fig. 2 for (3 = 0-20% and 
60-80%, that for 30-50% being in between the two. Clearly, the px distributions are affected 
by the increase of r\ mainly in the large-pr region. There are smaller changes at low pt as 
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Table 2: Parameters in Eqs. Q and (fTU|) 



y 


% 


9 


u 


d 


s 


u 


d 




A 


196.52 


55.65 


60.74 


3.114 


11.55 


12.23 


-0.75 


B 


1.442 


1.064 


1.045 


1.657 


1.330 


1.292 




n 


8.654 


7.533 


7.483 


8.798 


8.385 


8.319 




A 


254.06 


61.64 


65.26 


3.953 


13.35 


13.27 


-0.25 


B 


1.265 


0.996 


0.990 


1.456 


1.228 


1.223 




n 


8.207 


7.314 


7.293 


8.320 


8.102 


8.101 




A 


244.63 


59.51 


61.57 


3.786 


13.02 


12.50 


0.25 


B 


1.260 


0.991 


0.993 


1.453 


1.218 


1.229 




n 


8.175 


7.281 


7.280 


8.292 


8.068 


8.099 




A 


177.86 


51.35 


51.93 


2.745 


10.88 


10.34 


0.75 


B 


1.419 


1.039 


1.050 


1.646 


1.291 


1.308 




n 


8.546 


7.420 


7.439 


8.711 


8.263 


8.308 




A 


132.78 


43.64 


43.70 


2.030 


9.016 


8.580 


1.0 


B 


1.600 


1.103 


1.120 


1.874 


1.384 


1.402 




n 


8.959 


7.590 


7.621 


9.192 


8.496 


8.544 




A 


12460 


5.68e6 


1.77e7 


65.35 


6.562el2 


1.045e8 


2.2 


B 


0.3184 


(-)0.03396 


(-)0.0277 


0.4420 


0.2662e-3 


0.02054 


(k = 21.9) 


n 


5.939 


4.873 


4.897 


5.900 


5.191 


5.241 




m 


7.000 


5.320 


5.341 


7.657 


5.966 


5.951 




A 


10080 


2.391e5 


2.349e5 


30.396 


485.45 


814.80 


3.2 


B 


0.3360 


(-)0.05117 


(-)0.05218 


0.6214 


(-)0.1410 


(-)0.1287 


(ko = 8.13) 


n 


5.977 


4.539 


4.555 


6.545 


4.574 


4.607 




m 


6.024 


4.548 


4.559 


6.189 


4.943 


4.948 
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seen in the log scale, although they are not negligible in the linear scale, since the spectra 
there are proportional to C 2 ((3,r)) that varies substantially with rj according to Table I. 
The 7] dependence of the pt distributions in Fig. 2 is our prediction for which we have not 
adjusted any free parameters. No data on identified pions are currently available to check 
those results. 

The notable feature of Fig. 2 is that at rj = 3.2 the distributions behave nearly as 
straight lines in the log plot. The exponential behavior suggests that only the soft partons 
contribute to the pion formation. To see this more clearly, we show in Fig. 3 the different 
contributions to the spectra for f3 = 0-20% and rj = 3.2. Indeed, the TS and SS components 
are much smaller than the TT component, and are insignificant for px < 3 GeV/c. The 
reasons are twofold: not only is the soft parton density lower at 77 = 3.2, but also the hard 
parton distributions are severely suppressed at high kx- The former is evident in Table I; 
the latter is not as obvious in Table II. We plot the gluon distributions in Fig. 4 for the four 
values of 77, and see the precipitous fall for rj = 3.2, as k? approaches the kinematical limit 
at kx = 8.13 GeV/c. The suppression of fi(kT,i]) does not reduce the S term quadratically 
because fi(k T ,r]) appears only once in Eq. (J7|). Consequently, the TS and SS components 
can have comparable magnitudes in Fig. 3. Their significantly reduced contribution to the 
overall distribution exposes the TT contribution to be the dominant component forpr up to 3 
GeV/ c. One can reasonably question the validity of extrapolating the soft parton distribution 
T(pi) to pi ~ 1.5 GeV/c in its exponential form. Our view is that, instead of adopting some 
low-pr model that has its own ambiguities, it is sensible to use the exponential form of Eq. 
(j3J) for the soft parton distribution as a working hypothesis without introducing extra free 
parameters so as to make predictions that can be tested experimentally. The important 
observation is that the hard partons are suppressed at high 77 and that any prediction by 
pQCD should not neglect the soft background, which is shown to be more important than 
fragmentation at high 77. 

Having obtained the pion spectra at all j3 and 77, we can now calculate the central-to- 
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peripheral ratio 



Rcp(M 



dN n /pTdprdri(P)/ (iVcoiiOg)) 
dN n /p T dp T d v ({3 p )/(N coll ({3 p )) ' 



(11) 



where the reference (3 P = 60-80% and (iV co ii(/3)) is the average number of binary collisions at 
(3. The results for (3 = 0-20% and 30-50% are shown in Fig. 5 for the four values of 77. The 
data points are for (h + + h~)/2 for 77 = and 1, and for h~ for 77 = 2.2 and 3.2 [5]. The case 
of 77 = shows the Cronin effect that is well described by our result where the solid line is for 
(3 = 0-20% and the dashed line for j3 = 30-50%. As in Ref. [UJ, no k T broadening by multiple 
scattering in ISI has been put in. For 77 > 1, although the agreement of our results with 
the data is not perfect, they nevertheless exhibit the essence of the trend, i.e., Rcp becomes 
smaller as 77 is increased. That feature has been regarded as the distinctive characteristics of 
forward production, and is now approximately reproduced by our treatment that contains no 
new physics and no adjustable parameters. The case of 77 = 3.2 is the simplest to interpret, 
since the shower contribution is insignificant. The constancy of Rqp in our result for < 3 
GeV/c is a consequence of the fact that we have fixed the value of the inverse slope T for the 
soft parton distribution, independent of centrality. The suppression of Rcp at 77 = 3.2 is due 
to the decrease of C(/3,rj) with increasing rj and the insufficiently fast decrease of C(/3, 77) 
with increasing (3 to overcome the decrease of {N co n((3)) that rescales the spectra in Eq. ((TT} . 
In short, since the density of soft partons diminishes as one goes far into the deuteron side, 
less particles are produced by the recombination of those soft partons. 

In the foregoing we have fixed T for all (3 and 77 as an assumption for the sake of not making 
it an adjustable parameter. The theoretical results, as shown in Fig. 5, are remarkably close 
to the data. However, a mild dependence of T on (3 and 77 cannot be excluded. Indeed, the 
increase of Rcp with pt in the data at 77 = 3.2 suggests a decrease of T with (3. We adopt a 
simple parametrization of that dependence as follows 



T^, V )=T (l-e(3 V ). 



(12) 
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Since Rcp is plotted in linear scale in Fig. 5, it is possible to determine e by fitting the data 
at i] = 3.2, despite the absence of the spectra themselves (which would be in log scale and 
insufficiently accurate to determine small differences in T by themselves separately). We fix 
T = 0.208 GeV (more precisely T^ 1 = 4.8/GeV), and set (3 = 0.1,0.4, and 0.7 for 0-20%, 
30-50% and 60-80% centrality, respectively. We vary e to obtain the best fit of the data in 
open circles at i]=3.2 in Fig. 5, and get 

e = 0.0205. (13) 

The results for all other values of (3 and rj are shown in Fig. 6. Compared to Fig. 5, there is a 
slight improvement of the agreement with data at 77=1.0, but is a little higher than the data 
at 77=2.2 for p T > 2.5 GeV/c. Generally speaking, the trend of the data with increasing p T is 
better reproduced when T{f3, rj) is allowed to decrease slightly with (3 and 77. That decrease 
is less than 4.6% only even at the highest values of (3 (0.7) and 77 (3.2). Thus the constant T 
assumption is not a bad approximation and serves to reproduce the data reasonably well as 
in Fig. 5. However, we have not yet exhausted all aspects of physics that can influence the 
fit of the data. 

4 Momentum Degradation 

There is a piece of physics that we have not yet considered, but it is a phenomenological 
fact that should not be ignored. Baryon stopping generally refers to the loss of projectile 
proton momentum in pA collisions, as it passes through a target nucleus. Although such 
nomenclature is misleading from the point of view of the role that the proton constituents 
play, empirical evidence for the momentum degradation of the detected nucleon as a function 
of the nuclear size is not disputed. A number of experiments have shown that the produced 
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nucleoli distribution in pA collisions has the exponential form in xp 



dN N 



oc exp[— A(u) xp] , 



(14) 



dxp 



where the slope A(y) depends on the average number of collisions v |27| I28j. Baryon stopping 
loosely refers to the phenomenon that A(V) increases with v. Our question here is whether 
such a behavior has a dynamical origin that can affect our treatment of forward production 
in d+Au collisions. The question is relevant, since in both problems there is suppression of 
production probability at high 77. 

The production of leading nucleon in p + A collisions has been studied in the framework 
of the valon model for low-pr processes j2H|, and the distribution in Eq. (fTljl is obtained by 
attributing the momentum degradation effect to the projectile valons as they traverse the 
target nucleus [3UI- That is not contradictory to the information gained from the more recent 
experiments at RHIC, where energy loss effects are found to be absent at large pt in d+Au 
collisions pH} E21 Ell Ell • The former problem is at low pp and valons are dressed valence 
quarks of the proton, whereas the latter refers to hard partons that go through the nucleus 
at large angles with negligible interaction with the cold medium. It is the interpolation 
between these two extremes that is pertinent to the rj dependence in our problem here. 

The nucleon distribution in Eq. (fTl^l has been converted in the valon model to the pion 
distribution in the form 



where A = 0.2 |35j . There is no reliable way to relate that behavior in xp at low pt to the r\ 
dependence at intermediate pt without treating the transition from soft to hard processes. 
Since the boundary condition is that there is no energy loss at r\ = 0, we adopt the ansatz 



dN n 



oc exp[— \{u — 1) xp] , 



(15) 



dxp 
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that the degradation factor is 



C(A^) = exp[-«(JV c -l)7 ? ] , (16) 

where N c = (N co n((3)) and k is a parameter to be determined from A by matching Eqs. ([T5]) 
and ()16)1 at forward rapidity. The expression for (((3,rj) represents the property that the 
larger 7] is, the more time the constituents of the projectile spend in the valon state, while 
the valons propagate through the nuclear medium and suffer momentum degradation. 

To relate Eqs. ()15j) and (fTKJ) we note that at rj = 3.2, if (p?) = 2 GeV/c, the corresponding 
(xp) = 0.25. For (3 = 0-20%, we use the values v m 9 (for pAu collisions) and iV c « 15 (for 
dAu collisions), and get 

k«0.01. (17) 

We now use this value of k in Eq. (fTT)|) and apply (((3,7]) multiplicatively to the shower 
distribution S, but not to the soft distribution T , since the effect of degradation is already 
included in the determination of C((3,rj) through the use of the experimental values of 
dNch/ drj(P) in Eq. (JHJ). Equation (fT2"J) is used for T((3,rf). For SS recombination we do not 
apply (((3,7]) quadratically, since the shower partons are from one jet. 

With the degradation effect taken into account the results on the ratio Rcp are shown in 
Fig. 7. We note that there is improvement in the agreement with data compared to Fig. 6, 
especially at rj = 2.2. At rj = 1.0 the solid line no longer overshoots the dashed line at high 
Pt- There is no change at 77 = since (((3,0) = 1, and there is an improvement of the fit 
at 7] = 3.2. Since the data are for either (h + + h~)/2 or h~ , while our calculation is for 7r + 
specifically, perfect agreement between theory and experiment should not be expected. The 
effect of momentum degradation is at most 30% on the pr distributions at r\ = 2.2, so such 
changes are barely perceptible in the log plots of the spectra in Fig. 2, which therefore remain 
as the prediction of our treatment. Rcp in linear scale reveals the degradation effect more 
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sensitively. It is evident from Fig. 7 that the essence of forward production of pions in d+Au 
collisions is essentially captured in our description of hadronization by parton recombination, 
when T is allowed to depend on (3 and rj, and when the effects of momentum degradation is 
taken into account. 

5 Asymmetry Ratio for Backward to Forward Rapidi- 
ties 

So far we have restricted our study to only the forward region. The backward region on the 
gold side contains the properties of the nucleus not present at large 77, and should therefore 
behave differently from what we have obtained in previous sections. Recent data from 
STAR show significant asymmetry in the ratio of the charged hadron spectra for backward 
to forward rapidities in the range 0.5 < \rj\ < 1.0 [SHI- That ratio reaches a peak higher than 
1.3 for 2 < px < 3 GeV/c. At a qualitative level the phenomenon can easily be understood 
in our approach to the problem, since there are more soft partons at rj < than at rj > 0. 
We now want to examine the asymmetry quantitatively as another test of our treatment of 
hadronization. 

The formalism for particle production at rj < is the same as for rj > 0. Both T and 
S must, however, change, as rj enters the negative region. The parametrizations of the hard 
parton distributions fi(hr,y) are already given in Table II. Since the data for the backward 
region are for the range —1.0 < t] < —0.5, we shall use the parameters for y = —0.75 in 
Table II. Similarly, the values for y = 0.75 will be used for the forward region. For the soft 
parton distribution T(pi) we continue to use Eq. (JBJ to determine C(j3,rj) with dN c ^/ 'drj(P) 
taken directly from the data [23j. For f3 = 0-20%, we obtain C(J3, -0.75) = 12.372 and 
C(/3, 0.75) = 11.527 (GeV/c) -1 . This represents a small, but significant, asymmetry of the 
soft parton density. As for the 77 dependence of the inverse slope T(/3,r]) we use the formula, 
Eq. ()12j) . already determined in Sec. 3 from the region 77 > 0, now applied to rj = ±0.75. 
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We now can calculate the pr distributions for ir + using the appropriate values of C, T 
and fi at 77 = ±0.75. The results are, of course, not visually distinguishable from that at 
77 = in Fig. 2 plotted in log scale. However, their ratio plotted in linear scale is more 
sensitive to the small changes. The backward/forward asymmetry ratio is defined by 



dN v / 'dpT drj (77 = 


-M) 


dN n /dpr drj (77 = 


+\v\) 



Rb/f(pt, \v\) = ,, r 7~7~ -7-7— , i-fv • (18) 



The data for R B / F (p T , 0.75) as shown in Fig. 8 are for all charged hadrons Thus our 
calculation for 7r + alone is not enough for the purpose of comparison with the currently 
available data. In the same way that we have treated proton production in Au+Au [HI] 
and d+Au collisions |21| . we calculate the proton spectrum at 77 = ±0.75. Furthermore, we 
take the 7r~ yield to be the same as for tt + , and p yield to be 0.7 of that of p. The sum 
of 7r + + 7T _ ± p + p is shown by the solid line in Fig. 8. The result has the correct rise for 
Pt < 2 GeV/c, but is lower than the data at higher p T . We expect that the production 
of kaons can further increase the theoretical curve, but since they involve strange quarks 
that are enhanced in the soft component, we do not digress here to that peripheral subject. 
Our present result from the non-strange sector is sufficient to indicate that the asymmetry 
data can be understood in our approach to hadronization. Note that our result on Rb/f has 
been obtained without any new free parameter. A better way to compare theory with data 
would be to have identified pions at 77 = ±0.75, which is within the feasibility of some RHIC 
experiment. 

The conclusion that one can draw from this study of the backward-forward asymmetry 
is that there is no transition of basic physics from multiple scattering in ISI on the 77 < 
side to gluon saturation on the 77 > side [1. Our emphasis on the hadronization process 
in the final state provides a universal framework for the description of particle production 
at all 77, p, and pr- 
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Based on saturation physics, a recent calculation of the \ow-pr distribution of charged 
hadrons in minimum bias d+Au collisions at rj = 3.2 appears to have good agreement with 
data, when only the scattering of qq dipole on the nucleus described as CGC is taken into 
account [37] ■ However, the contribution of gluons to the cross section is not negligible, since 
the density of gluons in a proton at x — 0.25 is more than half the density of the quarks at 
Q < 5 GeV/c |38j . Furthermore, hadronization by use of fragmentation function is subject to 
the usual question about the p/ir ratio that has shown the inadequacy of the fragmentation 
model. Thus the agreement with data at low pt seems fortuitous. 

It has recently been stated that the Cronin effect seen at midrapidity goes away at 
forward rapidities because the supposed origin of the enhancement at rj — (i.e., multiple 
scattering in the initial state) is replaced by CGC that is responsible for the suppression at 
i] > [TJ El EZj - The results of this work combined with those of Refs. [TTJ EI] suggest that 
neither mechanisms are the primary causes of the enhancement and suppression, and that 
both phenomena are the consequences of the same hadronization process by recombination. 

6 Conclusion 

On the basis of parton recombination we have successfully described pion production at 
intermediate px in d+Au collisions in the forward rapidity region. The formalism is an 
extension of the one at midrapidity where the Cronin effect has been explained in terms of 
FSI only, and where the experimental fact R P CP > K^ P is interpreted as a consequence of 
the dominance of 3-quark recombination over fragmentation. In the extension to rj > only 
one new parameter, e, is introduced to describe the j3 and rj dependence of T, but no new 
physics has been added. The suppression of Rep at r\ > 1 is due mainly to the reduction of 
the density of soft partons that recombine either among themselves or with the semi-hard 
shower partons. The effect of momentum degradation responsible for baryon stopping has 
been considered, and is found to have a minor effect on Rcp, although it does render a 
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better agreement with the data at all rj and pr- The reduction of soft parton density in the 
forward direction may be related to gluon saturation, but there is no explicit reliance on 
small-x physics in the calculation. At large 77 the production of hard partons is suppressed, 
so pQCD calculations supplemented by fragmentation is likely to underestimate the pion 
spectra at intermediate pt- 

Extending the consideration to the backward region, we have used the same dependence of 
the inverse slope T((3, 77) on 77, now extrapolated from 77 > to rj < 0. We find that the general 
properties of the backward-to-forward ratio of the charged hadron distributions can be 
reproduced by our calculation of only the non-strange sector. Thus, particle production at 
all i] can be described by the same formalism for all pr at any centrality. There is no need 
for any change of physics in going from the backward to the forward region. 

To further verify the validity of our treatment, the predicted pion spectra in Fig. 2 should 
be checked by experiments. To have the correct Rcp, as in Fig. 7, is only the necessary 
condition for the underlying physics to be relevant, but not sufficient. Spectra themselves 
may disagree with the data, yet still have their ratio come out right. Proton spectra should 
also be measured at i] > and the p/n ratio shown as a function of pt, (3, and 77. Although 
we have not calculated the pt distribution of proton for 77 > 1, we expect that the p/n ratio 
would not be small for rj > 1, although lower than the maximum of the ratio for rj = 
because the lower soft parton density at 77 > 1 inhibits the formation of protons more than 
it does pions. The measurement of that ratio will provide a severe test on any model of 
particle production at any 77. 

Another area of investigation that can shed light on forward physics is to determine the 
presence or absence of back-to-back jets in azimuthal correlation when a particle is detected 
at large 77 and intermediate pr- If the pt distributions of the various components contributing 
to the pion spectrum at i] = 3.2 in Fig. 3 is correct, then we do not expect any significant 
jet signature until pr > 3 GeV/c. Even at pr = 4 GeV/c the dominant component is TT 
recombination, so at A0 = 71 there should be only a small jet-like component that stands 
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above a high level of background uniform in 0. Some aspect of that feature has already been 
observed in the preliminary data of STAR [3*H| . 

Since hadronization of partons by recombination is a process in the final stage of evolution 
of the partons, it is not in conflict with any dynamical model that correctly describes the 
beginning and subsequent evolution of those partons. Thus our work here provides the 
necessary link between the predicted parton spectra and the observed hadronic data for 
testing the validity of any proposed model. 
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Figure Captions 



Fig. 1. Transverse momentum distributions of ir + produced at midrapidity in dAu collisions 
for three centrality cuts. Light solid lines are for the recombination from the TT 
component, dashed lines for TS and dashed-dotted lines for SS components. Heavy 
solid lines are the sums of all three components. 

Fig. 2. Transverse momentum distributions of ir + produced at different pseudorapidities 
for two centrality cuts. 

Fig. 3. Transverse momentum distributions of 7r + produced at 77 = 3.2 for 0-20% centrality, 
showing the three components. 

Fig. 4. Distribution of gluons produced in central d+Au collisions at four pseudorapidities. 

Fig. 5. R CP for 0-20%/60-80% (filled circles and solid lines) and 30-50%/60-80% (open 
circles and dashed lines) for four pseudorapidities when T is assumed to be constant. 
Data are from No momentum degradation is considered in the calculation. 

Fig. 6. Same as for Fig. 5 but with T(j3,rj) given by Eqs. (fT^ and (fTT?|) . 

Fig. 7. Same as for Fig. 6 but with momentum degradation considered. 

Fig. 8. Ratio of p? distributions for backward to forward pseudorapidities at \r]\ = 0.75. 
Data are for all charged hadrons from 36J. Solid line is the calculated ratio for 7r + + 

7T~ + p + p. 
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